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I.  INTRODUCTION 


An  accurate  cold  beam  linear  dispersion  relation  for  the  negative  mass 
Instability  In  high  current  conventional  or  modified  betatrons  Is  critical 
for  designing  experiments,  benchmarking  computer  simulations,  and  develop¬ 
ing  warm  beam  dispersion  relations.  In  this  report  we  derive  a  dispersion 
relation  correct  to  first  order  in  the  toroidal  aspect  ratio  of  the  beta¬ 
tron  valid  for  small  toroidal  mode  numbers.  Limited  comparisons  with  ear¬ 
lier  work  are  provided  for  reference. 

The  two  earlier  treatments  of  negative  mass  Instability  growth  In  high 
current  betatrons  Included  toroidal  curvature  effects  In  the  particle 
(lynamics  but  not  In  the  electromagnetic  field  equations.^ In  other 
words,  the  dispersion  relations  were  derived  for  beam  motion  in  a  toroidal 
cavity  but  with  fields  computed  for  an  off-center  beam  in  a  straight  tube 
of  circular  crossectlon.  The  first  analysis,  that  of  Sprangle  and 
Vomvorldls,  used  an  equation  of  continuity  appropriate  to  a  straight  cyl¬ 
inder,  as  well.  In  contrast,  the  second  analysis,  due  to  Hughes  and 
Godfrey,  Incorporated  an  equation  of  continuity  correct  for  a  toroidal 
beam.  Predicted  growth  rates  from  these  two  models  often  have  been  found 
to  disagree  by  factors  of  two  for  parameters  corresponding  to  peak  growth 
and  by  even  more  at  higher  beam  energies.  More  disturbingly,  both  models 
systematically  overestimate  Instability  growth  rates  observed  in  three- 
dimensional  computer  simulations.^ ***  We  attribute  these  discrepancies 
to  Incomplete  treatment  of  toroidal  curvature  effects. 

The  dominant  curvature  correction  to  the  field  equations  is  easily 
Identified.  Radial  motion  of  an  electron  ring  (see  Fig.  1)  gives  rise  to 
azimuthal  bunching.  This  constitutes  the  difference  between  charge  contin¬ 
uity  In  cylindrical  and  toroidal  geometry.  Conversely,  azimuthal  bunching 
gives  rise  to  a  radial  electric  field,  which  constitutes  the  (main)  differ¬ 
ence  between  electromagnetic  fields  In  cylindrical  and  toroidal  geometry. 


Schematic  representation  of  betatron  cavity  with  electron  beam 
Mirror  and  optional  toroidal  magnetic  fields  are  applied. 


The  additional  coupling  between  radial  and  azimuthal  beam  motion  need  not 
be  small.  Indeed,  for  sufficiently  large  energies  this  coupling  through 
the  fields  must  dominate  coupling  through  the  particle  dynamics.  The 
latter  falls  off  with  energy,  while  the  former  does  not. 

The  present  paper  includes  these  and  other  (lesser)  curvature  effects 
to  first  order  in  the  equilibrium  and  perturbed  fields,  producing  what  we 
believe  to  be  a  more  accurate  dispersion  relation.  Like  its  predecessors, 
this  analysis  ignores  the  internal  dynamics  of  the  beam,  treating  it  as  a 
string  of  rigid  disks.  The  assumption  of  long  azimuthal  wavelengths  also 
is  implicit  in  our  derivation.  Relaxing  the  last  constraint  is  conceptu¬ 
ally  straightforward  but  algebraically  difficult. 

In  Section  II  we  assemble  the  various  particle  and  field  equations 
into  a  5x5  self-adjoint  matrix  operator  acting  on  the  perturbed  beam  cen¬ 
troid  position  and  the  electrostatic  and  electromagnetic  potentials.  The 
differential  equations  are  solved  to  first  order  in  the  toroidal  aspect 
ratio  in  Section  III,  and  the  remaining  algebraic  equations  collapsed  to 
the  desired  dispersion  relation.  When  no  toroidal  magnetic  field  exists 
(i.e.,  for  a  conventional  betatron),  an  analytic  growth  rate  expression  for 
the  negative  mass  instability,  valid  over  a  wide  parameter  range,  can  be 
derived  from  the  general  dispersion  relation,  and  this  is  done  in  Section 
IV.  We  have  not  yet  attempted  to  obtain  a  corresponding  simple  growth  rate 
expression  for  the  modified  betatron. 

Finally,  Section  V  presents  a  reevaluation  of  negative  mass  growth  for 
the  Office  of  Naval  Research  racetrack  induction  accelerator  design,^  based 
on  the  new  dispersion  relation.  The  average  instability  growth  rate  during 
acceleration  from  a  few  to  50  MeV  is  reduced  by  a  factor  of  two  relative  to 
an  earlier  prediction.^  With  this  improvement,  the  proposed  accelerator 
may  be  able  to  "outrun"  the  instability  with  only  a  modest  beam  thermal 
energy  spread.^ 

Comparisons  with  computer  simulation  results,  as  applied  to  the  Naval 
Research  Laboratory  modified  betatron  design,  are  described  elsewhere.^ 
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II.  FORMULATION  OF  THE  PROBLEM 


We  consider  an  electron  ring  circulating  in  an  azimuthally  symmetric 
conducting  cavity.  Dimensions  and  coordinates  are  depicted  in  Figure  1. 

For  the  present  we  assume  only  that  the  beam  and  cavity  crossections  are 
symmetric  about  z=0.  Otherwise,  the  crossections  may  assume  any  reasonable 
shape.  The  characteristic  minor  radius  of  the  beam  must  be  much  less  than 
that  of  the  cavity  so  that  the  internal  dynamics  of  the  beam  can  be 
ignored.  A  mirror  magnetic  field  maintains  the  beam  in  a  circular  orbit, 
and  a  toroidal  field  can  be  added  to  confine  the  beam  against  its  self¬ 
fields.  which  need  not  be  small. 

Linearized  equations  for  small  transverse  displacements  of  the  beam 
centroid  are  easily  obtained  from  the  single  particle  equations  of  motion. 


Perturbed  quantities  are  preceded  by  a  delta  (e.g.,  6z),  while  unperturbed 
quantities  are  not.  Total  time  derivatives  of  perturbed  particle  quanti¬ 
ties  are  represented  by  dots  above  the  quantities  (e.g.  6z).  Note  that  the 
last  term  in  Eq.  (2)  can  be  rearranged,  if  desired,  by  means  of  the 
equilibrium  radial  force  balance  equation. 


Although  the  perturbed  azimuthal  velocity  equation  can  be  derived  In 
many  ways,  using  the  single  particle  energy  formula  seems  simplest. 
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The  perturbed  azimuthal  angle  <58  of  a  beam  disk  Is  related  to  <5Ve  by 
6Vq  =  R  66  +  6r  (5) 

Perturbed  beam  currents  resulting  from  the  centroid  displacements 
described  by  Eq.  (1),  (2),  (4),  and  (5)  are 

6J^  =  p  6r,  6J^  =  p  6z 

and 

=  P  fiVg  +  Vg  6p  (6) 

Perturbed  charge  density  Is,  In  turn,  derived  by  substituting  the  perturbed 
currents  Into  the  continuity  equation  and  integrating  the  result  in  time. 

6p  +  1^  p60  +  1  Ip  r  p6r  +  p6z  =  0  (7) 

Alternatively,  Eq.  (7)  can  be  obtained  by  considering  how  the  density  of  a 
beam  element  changes  as  It  Is  displaced  Infinitesimally  In  each  direction. 

As  explained  In  the  Introduction,  we  are  limiting  consideration  to  low 
frequencies  and  long  azimuthal  wavelengths.  In  this  limit  we  need  only 
determine  the  electrostatic  and  azimuthal  electromagnetic  potentials,  which 


(9) 


i  Ip  r  +  6A  =  -  6Jg 

The  electric  and  magnetic  field  combinations  appearing  in  Eq.  (1),  (2),  and 
(4)  are  expressed  in  terms  of  these  potentials  as 


6E^  -  Vg  =  - 

It 
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(10) 

«Ep  +  V3  6B^  =  - 
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These  equations  complete  our  model. 
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Let  us  now  cast  the  equations  in  matrix  form  for  compactness  and  to 
emphasize  their  symmetry.  In  so  doing  we  also  Fourier  transform  the  equa¬ 
tions  in  time  and  azimuthal  angle,  i.e., 

n  *  iQ  * 

and  eliminate  6V0  and  the  perturbed  electric  and  magnetic  fields. 
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Note  that  the  matrix  operator  Is  self-adjoint,  as  one  would  hope.  Certain 
symbols  appearing  in  Eq.  (13)  are  defined  as 


Equation  (13)  can  be  solved  numerically  for  beams  and  cavities  of 
arbitrary  cross  section.  (Before  doing  so.  It  probably  would  be  desirable 
to  perform  a  conformal  transformation  on  r,2  to  map  the  beam  and  cavity 
boundaries  onto  coordinate  surfaces.)  In  the  next  Section  we  Instead 
assume  the  beam  and  cavity  to  be  circular  and  concentric  in  cross  section 
as  Illustrated  In  Fig.  1,  In  order  to  obtain  an  analytical  disc  .’on 
relation.  We  should  bear  In  mind,  however,  that  other  configi  tions  may 
exhibit  improved  stability. 


III.  DERIVATION  OF  THE  DISPERSION  RELATION 


To  develop  a  dispersion  relation  from  Eq.  (13),  even  for  the  simple 
geometry  explicitly  represented  by  Fig.  1,  is  a  long  and  involved  process. 
For  the  sake  of  brevity,  we  here  restrict  ourselves  to  outlining  the  pro¬ 
cedure  and  citing  pertinent  results. 

Converting  the  partial  differential  operators  in  Eq.  (13)  to  a  set  of 
coupled  ordinary  differential  operators  is  the  first  step.  To  do  this  we 
carry  out  the  coordinate  transformation 

z  =  X  cost|/,  r  =  R  +  X  Sint 

and  expand  the  potentials  as  Fourier  series  in  poloidal  angle 

fit  =  I'5<|i"'(x)  e’’"'*',  6A  =  I6A"'(x) 

Next,  we  recast  Eq.  (13)  as  a  variational  integral,  insert  the  6(t>  and  fiA 
expansions,  explicitly  integrate  over  t.  and  perform  variations  with 
respect  to  and  6A"'.  This  results  in  an  infinite  set  of  different¬ 
ial  equations  in  ^  coupling  all  the  poloidal  modes  of  the  potentials. 

Expanding  the  differential  equations  in  x  with  respect  to  the  cavity 
aspect  ratio,  b/R,  leads  to  a  natural  truncation  of  the  infinite  system. 
Basically,  an  expansion  to  order  |m|  in  b/R  is  consistent  with  dropping  all 
higher  poloidal  mode  numbers.  Two  orderings  spring  to  mind.  In  the  first 
0),  i/R  and  the  electron  oscillation  frequencies  all  are  taken  to  be  of 
order  unity  (or  smaller),  and  the  equations  expanded  to  order  b/R.  In  the 
second,  to  and  i/R  are  taken  to  be  of  order  b/R,  and  the  equations  expanded 
to  b^/R^.  The  first  alternative  is  much  simpler,  seems  to  capture  the 
essential  physics,  and  agrees  reasonably  well  with  simulation  results,  so 
we  use  it  in  the  following  calculations.  We  are,  however,  investigating 
the  consequences  of  the  second  ordering  and  will  publish  our  findings  at  a 
later  date. 


With  the  ordering  selected,  straightforward  but  tedious  calculations 
yield  for  x  <  a: 
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The  validity  of  these  solutions  to  the  truncated  equations  has  been  veri¬ 
fied  using  the  symbolic  manipulation  program  MAXIMA. 

Back  substitution  of  the  potentials  into  the  variational  integral 
allows  the  remaining  integrations  to  be  carried  out  explicitly.  (For  con¬ 
sistency,  the  results  are  truncated  to  first  order  in  b/R.)  Again,  MAXIMA 
was  used  to  check  the  manipulations.  The  resulting  matrix  equations  for 
the  perturbed  centroid  location  appear  as  Eq.  (20).  Its  determinant  is 
the  sought-for  dispersion  relation. 
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In  evaluating  the  matrix  elements,  it  is  necessary  to  know  the  equi¬ 
librium  fields.  They  are^ 

Ep-  7  ('•-'»*  TO  (E‘) 

*  Bz°  (‘  -  "  (22) 

with  similar  expressions  for  £2  and  Bp.  Eq.  (3)  determines  the  magni¬ 
tude  of  Bz®,  and  Be®  is  arbitrary. 

The  dispersion  relation  is 

b“^/y^  =  0  (23) 

It  represents  two  longitudinal  (m  =  0)  modes,  described  by  the  longitudinal 
dielectric  constant, 

^  ^  (i  ^  2  -  “^)  (24) 

and  four  transverse  (m  =  ±1)  modes,  described  by  Eq.  (23)  with  x  =  0* 
transverse  oscillation  frequencies  uz  and  up  are. 
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Ep  and  Bz  in  Eq.  (26)  are  given  by  Eq.  (21)  and  (22)  with  r=R  and 
Bz®  omitted.  The  last  term  in  Eq.  (26)  is  very  small.  Typically,  the 
external  field  index  n  is  chosen  such  that  Note  that  v  is 

Budker's  parameter,  equal  to  pa^/4. 

The  key  result  of  our  analysis  is  the  coupling  coefficient  between 
longitudinal  and  transverse  modes,  which  determines  the  negative  mass 
instability  growth  rate. 
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For  comparison,  the  dispersion  relations  from  Ref.  1  and  2  also  can  be  cast 
in  the  form  of  Eq.  (23)  but  with  coupling  coefficients,  respectively, 


and 


X  = 
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It  should  be  noted  that  there  is  a  degree  of  arbitrariness  in  the 
choice  for  the  functional  form  of  p(r,z).  In  the  preceding  analysis  we 
assumed  p  to  be  constant  out  to  the  beam  minor  radius  a,  where  it  drops 
abruptly  to  zero.  One  might  instead  have  chosen  pr/R  to  be  constant,  for 
instance.^  Fortunately,  such  changes  lead  only  to  insignificant  modifica¬ 
tions  of 


IV.  Be®  =  0  GROWTH  RATE  FORMULA 


In  the  absence  of  a  toroidal  magnetic  field,  the  negative  mass  insta¬ 
bility  dispersion  relation  reduces  to 

-  x/e  =  0  (30) 

A  simple  analytical  growth  rate  expression  can  be  derived  from  Eq.  (30)  in 
a  reasonably  straightforward  manner,  if  is  much  less  than  and  can 
be  dropped  from  the  equation.  We  address  the  validity  of  this  assumption 
at  the  end  of  this  Section. 

Solely  for  the  sake  of  algebraic  simplicity,  we  invoke  four  approxima 
tions. 


All  four  are  well  satisfied  for  typical  cases  of  interest.  Then, 


significant  cancellations  occur  when  Eq.  (32)  Is  expanded. 


(33) 


The  term  linear  In  £l.  In  fact,  only  survives  due  to  the  difference  In  the 
factors 


1  +  2An  I  ,  I  +  2An  j 

which  appear  In  the  first  and  second  terms,  respectively,  of  the  definition 
for  X*  Modifications  to  the  expression  for  x  discussed  near  the  end  of  the 
preceding  Section,  due  to  changing  the  functional  form  of  p,  are  of  higher 
order  than  the  terms  In  Eq.  (33)  and  so  drop  out.  Nonetheless,  we  should 
not  be  surprised  If  more  accurate  dispersion  relation  derivations,  perhaps 
Including  high  frequency  electromagnetic  effects,  adjust  the  linear  term 
somewhat. 

Substitution  of  Eq.  (31)  and  (33)  Into 

2  ^  n 

(o„  e  +  X  “  0 
r  ^ 


leads  Immediately  to  the  desired  growth  rate  formula. 
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In  obtaining  Eq.  (34),  we  have  taken 


V  2 

■  (!-")  -4  . 

which  Is  accurate  for  not  too  small,  an  assumption  we  have  already 
made.  Equation  (34)  reproduces  the  exact  numerical  solutions  of  the  dis¬ 
persion  relation  presented  In  the  next  Section  to  an  accuracy  of  10%. 

Equation  (34)  predicts  an  Instability  cutoff  for 

vy  >  4(1  -  n)  +  2An  (35) 

At  smaller  values  of  vy,  the  negative  mass  Instability  reduces  to  the  well 
known  expression^® 

2 /2 

^  *  2“  j)]  (36) 

We  remark  for  completeness  that  Eq.  (28),  from  Ref.  1,  and  Eq.  (29),  from 
Ref.  2,  also  yield  high  vy  cutoffs.  In  the  case  of  Eq.  (29)  Instability 
ceases  for^’^^ 

vy  >  4(1  -  n)/(l  +  2in  j)  (37) 

while  Eq.  (28)  has  the  same  limit  without  the  factor  of  four.  Both  analy¬ 
ses  lead  to  Eq.  (36)  for  sufficiently  small  vy. 

Let  us  return  to  our  Initial  assumption. 


V.  THE  ONR  RACETRACK  INDUCTION  ACCELERATOR 


The  Office  of  Naval  Research  together  with  the  Naval  Research  Labora¬ 
tory  have  developed  the  design  for  a  racetrack  induction  accelerator 
intended  to  accelerate  a  1  kA  (v  =  0.0588)  electron  beam  to  40  MeV  (y  =  80) 
in  fifty  revolutions.  In  an  earlier  report,^  we  investigated  the  beam 
breakup  and  negative  mass  instabilities  for  this  design,  concluding  that 
the  beam  breakup  instability  poses  few  problems  for  induction  module  O's  of 
order  six,  a  realistic  value.  The  negative  mass  instability  was  found  to 
be  somewhat  more  threatening,  although  there  was  reason  to  hope  that 
thermal  effects  would  reduce  the  predicted  growth  rates,  especially  for 
large  toroidal  mode  numebrs.^  Our  work  was  based  on  the  dispersion  rela¬ 
tion  of  Ref.  2.  Here,  we  repeat  the  analysis  using  Eqs.  (23)  -  (27). 

We  treat  the  racetrack  cavity  as  a  torus  of  major  radius  R  =  70  cm. 
(Straight  sections  of  the  racetrack  are  expected  to  have  a  favorable,  but 
very  small,  effect  on  total  negative  mass  instability  growth.)  Varying  the 
assumed  radius  has  little  affect,  because  changes  in  the  growth  rate  are 
approximately  balanced  by  changes  in  the  path  length.  In  accordance  with 
the  accelerator  design,  the  cavity  and  beam  minor  radii  are  taken  to  be  b  = 
7  cm  and  a  =  1  cm.  A  Be"  =  2  kG  (wc  *  1.73)  guide  field  is  applied. 

The  numerically  determined  growth  rate  for  these  parameters  and  l  = 

13,  a  toroidal  mode  number  near  the  upper  end  of  the  range  for  which  the 
model  is  valid,  is  plotted  as  the  solid  curve  in  Fig.  2.  The  result  from 
Ref.  6  is  shown  as  a  dashed  line.  The  average  growth  rate  is  reduced  by 
about  a  factor  of  two.  Our  more  optimistic  findings  suggest  a  total  nega¬ 
tive  mass  instability  amplification  factor  of  fifteen  e-foldings.  This 
growth  level  probably  can  be  reduced  to  an  acceptable  level  (say,  five  e- 
foldings)  by  thermal  effects  or  design  changes. 
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Figure  3  illustrates  the  corresponding  89®  =  0  growth  rates.  Peak 
values  are  reduced  somewhat  relative  to  the  earlier  results,  but  the  high 
energy  cutoff  is  shifted  up  in  energy  by  a  factor  of  twenty.  For  y  :>  20, 
Eq.  (34)  reproduces  the  solid  curve  well. 
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